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Abstract 



The rates of inclusive photo- and electroproduction of the 77c meson: 7(*) + p ^ rjc + X are 
calculated in the triple Regge region, integrated over the diffractive mass X. For the Regge 
] exchanges we use the hard pomeron and odderon, both being calculated in the framework 

■ of perturbative QCD. The integrated cross-section depends upon the coupling of the BFKL 

, pomeron to two C — —1 odderons, and it is found to be of the order 60 pb for photoproduction 

and 1.5 pb at = 25 GeV^. 

^ ; 1 Introduction 

The existence of the odderon [1], the partner of the Pomeron which is odd under charge conjugation 
^ ■ C, is an important prediction of perturbative QCD. It is a direct consequence of the number of 

colours Nc being greater than two. In the leading order, the odderon appears as a bound state of 
three reggeized gluons. Its experimental observation is a strong challenge for the experimentalists. 
A particular promising scattering process where the exchange of the odderon may be seen is the 
diffractive production of particles with a C-odd exchange, such as photo- and electroproduction of 
pseudoscalar mesons (PS), provided a large momentum scale is involved, which gives a justification 
for the use of perturbative QCD. This includes, in particular, the diffractive production of charmed 
pseudoscalar mesons, for example the process 7 + p — > ??c + Correspondingly, a large amount 
of literature has been devoted to this class of diffractive processes. For large photon virtualities 
Q^, for heavy mass PS mesons (such as rjc), and for large momentum transfers the relevant impact 
factors for the transition 7(7*) — >PS have been be calculated perturbatively [2]. As to the odderon 
structure, different models have been used: the exchange of three noninteracting gluons in a C=-l 
state [2], and, more recently, the perturbative QCD odderon with intercept exactly unity [3, 4], 
used to calculate the production rates of the rjc in [5] ^. In both models, there is some uncertainty 
coming from the coupling of the odderon to the proton. Numerical estimates for the cross sections 
turn out to be somewhat different in these approaches. However in all cases they are very small 

different approach, a nonperturbative odderon based upon the idea of a "stochastic QCD vacuum" has been 
used in [6]. 
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and, most unfortunately, do not grow with energy (in the case of the pcrturbative QCD odderon, 
they even slowly decrease with energy). This leaves little hope to see the odderon by raising the 
energy of the reaction. 

However the situation may become different if, instead of the quasielastic process ■~i+p f]c+P, 
one considers the inclusive cross section j + p ^ ijc + X in the triple Regge region. In this 
case the odderon does not couple directly to the quarks of the target proton but rather to the 
diffractive system 'X' which, for high masses, can modelled by a cut gluon ladder, the gluon 
density inside the proton. The proton is therefore coupled to the cut gluon ladder, i.e. the Pomeron, 
and this coupling is known through the gluon density. This fact permits to avoid the previously 
mentioned uncertainties in the odderon-proton coupling. Together with this process also the low 
mass diffractive state (the proton) with the meaning of a double odderon exchange is usually 
considered. 

In the Regge language this new situation basically involves the coupling of two odderons to a 
cut Pomeron, the POO vertex. Since we are using perturbative QCD both for the odderon and for 
the cut Pomeron, also this vertex has to be calculated in perturbative QCD. This has been done 
in [8]: the vertex has been obtained in an analysis of a six gluon amplitude Dq. In our application 
of this vertex we shall restrict ourselves to the leading large-Nc limit, which leads to a relatively 
simple form of Dq. In [8] it was also shown that the full amplitude Dq can be decomposed into the 
sum of two contributions, where the first one results from the reggeization of the gluon and , and 
the second one contains the POO vertex. Correspondingly, also our cross section comes as the sum 
of two pieces (denoted by P and POO, respectively). The second one corresponds to the normal 
'triple Regge picture' where the Pomeron splits into two odderons, whereas the first one is related 
to reggeization of the gluon and leads the exchange of three noninteraction gluons in the odderon 
channel. We will calculate both of them. For the odderon states we will use the solution found in 
[3], which has a maximal intercept (unity) and a very simple analytical form. This solution has 
already been used by us to calculate the odderon exchange in the process 7 +J9 — > ?7c + P [5]. 

The use of these elements allows to compute the diffractive cross section ^^l^i for the process 
7 + p — > r/c + -'^^ Our perturbative QCD analysis only depends upon one free parameter, the 
coupling of the gluon ladder to the proton: this coupling will be fixed by fitting the model to the 
gluon density of the proton. To simplify the calculations we restrict ourselves to the integrated 
(over M^) cross section: the calculation of the differential (in M^) cross section requires a slightly 
different treatment of the Dq amplitude which will not be pursued in the present work. Nevertheless 
the most inportant and basic information is given by the integrated cross section. 

As we have indicated before, we expect that the cross section for the inclusive process 7 + p ^ 
r\c^ X is larger than that for the quasielastic process 7 + p ^ r/c + P- We know that the cut 
gluon ladder grows as exp Ay, where y is the rapidity gap between the proton and the POO vertex, 
and A is the value of the Pomeron intercept minus unity. Prom this it follows that the bulk of 
the inclusive cross-section will come from the region where y is as large as possible, i.e. close to 
the total rapidity of the process (note however, that, in order to see the exchange of an odderon, 
one needs also a large rapidity gap between the outgoing r\c and the diffractive system). In other 
words, the mass of the diffractive system "X" wants to become as large as possible. Because of 
this growth (with energy) of the cut gluon ladder we expect to see a strong enhancement of the 
inclusive cross-section at high energies, compared to the quasielastic process j + p ^ r]c + p where 
the gluon ladder is absent. The comparison of our results with the quasielastic cross section is 
made difficult by the intrinsic uncertainty of the odderon-proton coupling: a recent analysis shows 
that the estimate obtained in [2] and also adopted in [5] may have used a too large value of this 
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coupling and has to be reduced: if this is the case, the cross section obtained in the present paper 
is, in fact, much larger that the quasiclastic one. 

Our paper is organized as follows. In the next section we shall briefly recall some results which 
constitute our starting point to attack the problem, the cross section formula. In section 3 the first 
contribution (P) is considered, and the corresponding integrated (in M^) cross section is written in 
terms of a multidimensional integral which, later on, will be be computed numerically. In sections 
4 and 5 the structure of the second contribution (POO) is considered. Both contributions are 
calculated in the large Nc limit, which leads to considerable simplifications and shows a symmetry 
shared by the leading odderon states [3] . Finally, the numerical analysis is presented and discussed 
in section 6, followed by the conclusions. 

2 The cross-section formula in perturbative QCD 

We start from the analysis [8] of QCD Feynman diagrams in the leading log s approximation, 
and we recapitulate the main results. The approach taken from [8] is a generalization of a previous 
analysis [9] of the 4 gluon system (related to the triple pomeron vertex): it extends this analysis 
up to 6 gluons in the t channel, and so it encounters, for the first time, the two-odderon state. As 
described in [8], the analysis of Feynman diagrams in the high energy limit leads to a tower of 
gluon amplitudes, D2, -D3, -D4, -D5, and Dq, which satisfy a set of coupled integral equations. These 
functions are non-amputated, i.e. they contain reggeon denominators for the outgoing (reggeized) 
gluon states (see Ref. [9]). The latter are more convenient degrees of freedom than the elementary 
gluons in this kinematics. In the present context we are interested in the Dq amplitude which can 
be used to build the cross section, integrated in the diffractive mass. 

In the analysis in refs. [8] and [9], all functions -D2, ••• start from the impact factor of a virtual 
photon which splits into a quark-antiquark pair. In the present case, the external particle is the 
proton: assuming that, to a good approximation, the proton can be viewed as a quark-diquark 
system, the coupling of the gluons to the proton should have the same structure as in the photon 
case; only the overall normalization of this coupling has to be treated as a phenomenological 
parameter. 

The diagrammatic structure of the differential cross section ^^^^j2 for the process 7+jo ^ rjc + X 
is illustrated in Fig. la. In the two exchange channels one recognizes the two odderon states, 
consisting of three gluons with pairwise interactions. The structure of the blob (related to Dq) 
will be discussed in a future paper. When the integration over the squared missing mass is 
performed, the expression for the cross section ^ = / dM"^ dtdM^ simplifies. The result is illustrated 
in Fig. lb: the blob now stands for Dq which can directly be taken from [8]. It depends upon the 
angular momentum variable ui which is conjugate to the total rapidity Y. 

The cross section can be written as 

t =^ E / ^-""V ^|i^^^(l,2,3)[$^(4,5,6)ri)6(l,2,3,4,5,6;u;). (1) 

Here denotes the impact factor for the transitions 7* r]c with the photon polarizations i = 1,2 
and the color structure date- Y is the overall rapidity; t = —q^ is the invariant associated to the 
momentum transfer across the impact factor, and the arguments 1,2,3 and 4,5,6 refer to both the 
colour indices Oj and transverse momenta, ki, of the gluons exchanged in the initial amplitude 
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Figure 1: Illustration of the process 7* + j3 — > ?7c + X. Diffractive cross section differential (a) and 
integrated (b) in the diffractive mass. 



i = 1,2, 3 and final (conjugated) amplitude i = 4,5,6. Finally, dfii is the integration measure for 
the 3 t-channel gluons on the Ihs: 

dm = d^kid^k2d^k^5'^{ki +k2 + k^-q) (2) 

and d/X2 is the analogous integration measure for the three t-channel gluons 4,5,6 on the rhs. The 
normalization factor ^ will be discussed in the next section. 

From the analysis [8] of the coupled integral equations it follows that can be presented 
as a sum of different terms. One term (denoted by D^) is obtained by collecting the rcggeizing 
pieces: the outgoing six gluon state may contain configurations where a pair of two gluons is 
in an antisymmetric color octet configuration, which satisfies the BFKL bootstrap condition and 
collapses into a single gluon. As a result, one obtains contributions with a smaller number of 
reggeized gluons. It is convenient to separate these configurations from the rest, i.e. to define the 
sub-amplitude which is 'irreducible' with respect to this reduction procedure. This reduction 
leads to the decomposition Dq = + separating the reggeizing (R) and irreducible (I) parts. 
The Dg term (cq.(6.3) of [8]) is rather lengthy; however, for an odderon in the (123) and (456) 
channels, we will need, in the large- A^c limit, only one term, denoted by W , which describes the 
transition of two reggeized gluons into six reggeized gluons: all other terms will be shown (see 
section 4) to be suppressed by a factor of 1/iV^ (or even higher powers of this). It is this piece of 

which yields the POO vertex. 

Diagrammatically, the piece of Dq which, in the large- Ag limit contains the Pomeron odderon 
vertex has the structure shown in Fig. 2. The internal blob - with two gluons entering from above 
and six gluons leaving below - defines the Pomeron odderon (POO) vertex, and its color structure 
is quite simple: 

'5fe,b'daia2a344aBa6^(l,2,3|4,5,6) , (3) 

where the 6, b' are the color labels of the reggeized gluons of the ladder above the POO vertex, 
Oj the color indices of the reggeized gluons below the vertex (counting from left to right). The 
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Figure 2: Illustration of the process ^* + p ^ r]c + X: the blob in the center of the figure denotes 
the transition vertex Pomeron — ^ 2 odderons. 



arguments of the function W refer to the momenta of the gluons. Below the POO vertex, we have 
the two noninteracting odderons: the pairwise interactions inside (123) and (456) lead to the color 
singlet odderon Green's functions. We note that this simple form emerges only after taking the 
large-iVc limit. In the more general case of finite Nc, the expression (3) has to be summed over 
permutations of the indices (123456). Moreover, in Fig. 2 below the POO vertex, we would have 
to include all pairwise interactions between the reggeized gluons. It is only in the large- A^'c limit 
that any rung which connects the two color singlet (123) and (456) costs a suppression factor of 
the order and, therefore, can be neglected. 

The Dq term is nothing but a sum of BFKL ladders in which, at the lower end, the reggeized 
gluons split into two three or four elementary gluons. Inserting this sum into the blob in Fig. lb 
and taking the large- AT^ limit, we arrive at structures illustrated in Fig.3: the BFKL ladder couples 
to odderon states consisting of three noninteracting gluons. Below we will discuss this in further 
detail: starting from the color structure of Dq, given in [8], eq.(6.2), it can be shown that all these 
contributions are subleading in l/Nc, in our further discussion we will keep one of them which is 
suppressed by a factor 1/Nc (all others are further suppressed). 

All this discussion refers to Fig. lb which illustrates the integrated inclusive cross section. In 
order to derive the pQCD formula for the differential cross section, an alternative separation of Dq 
is more suitable. For the case of D4 which leads to the triple Pomeron vertex, such a separation has 
been discussed in [9] and [10]: as a result, a slightly different expression for the triple Pomeron 
emerges. This question will be discussed in a forthcoming paper. 



Figure 3: The second contribution to the same process as Fig.l: the Pomeron couples to the two 
odderons which consist of three noninter acting gluons. 



3 Contribution of the reggeizing piece Dq 

Let us now analyze the two contribution to our inclusive cross section in some detail. We start with 
the ideologically (but not calculationally!) simpler part of the transition rate corresponding to the 
reggeizmg piece, D^, of [8] (Fig. 3). In the following, we will denote this piece by the superscript 
P. The corresponding inclusive cross-section is given by the expression 

^ = e E / ^-"V ^$^(l,2,3)[$^(4,5,6)ri;|^(l,2,3,4,5,6;a;). (4) 

The function Dq depending on all gluonic momenta is the reggeizing piece of the 6-gluon amplitude 
found in [8]. It is given by the sum of BFKL pomerons depending on various partial sums of the 
momenta of the 6 gluons 1,..6, multiplied by certain colour factors. All colour factors are obtained 
by permutations of inital (123) or final (456) gluons from 

^aiaiasaiasae ^ ^j. (^01^02^03^04^05^06) + ^ j. (^06 ^05 ^04 ^03 ^02 ^01 ) ^ 

where t"" is the quark colour matrix. This evidently gives 8 different colour factors, so that 
contains 8 terms with different colour structure. Coupling to the two impact factors in (1) 
one has to contract its colour indeces with a product (^010203 (^040506 corresponding to the C = — 1 
exchange of three gluons. Then one finds that all colour factors in are transformed into the 
same common colour factor Fc^^ 



J010203040506 J J 

(iV2-l)2(iV2_4)2 _ 

"'Oi0203"'a4a5a6 ~ SN"^ — c-^c ■ 
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At large Nc it is in correspondence with the general rules of the 1/Nc expansion. Note 

however at A'^c = 3 its value 200/27~ 7.4 is nearly 4 times smaller than given by the Nc — > oo limit 
243/8 ~ 30. Finally the overall factor ^ is equal to 



2 167r3 V4-4(27r)6(3!) 

The first factor 1/2 corresponds to the averaging over the two photon polarizations (we just consider 
transverse photon). Then one has the standard l/(167r^) phase space volume for the diffractive 
process times a tt factor due to angular averaging. In the squared term there is the contribution of 
(27r)~^ absent in each d/j, and an extra (27r)^^ which we associate to the impact factor [2] in our 
normalization. Moreover one has the symmetry factor 1/3! for each of the gluon triplets, 1/4 from 
the colour factor which is in fact {l/4:)dabc and 1/4 from the integrations over Sj, instead of over 
ki- in the definition of the impact factor. 

Separating the common factor g where gs is the strong coupling constant, we shall 

reproduce the rest of the amplitude from [8] in a simplified manner, taking into account that, 
first, in the pomeron of Fig. 3 the total momentum of the two gluons is zero and, second, in the 
high-energy limit the amplitude for the leading contribution is symmetric in the two gluons. Thus 
this ampitude P{k) (amputated, that is, without external gluon propagators) depends only on one 
of the gluon momenta. In terms of P{k) the reggeizing piece is then given by a sum of 31 terms: 

D^ = j2m- E pm- E ^'(i^)-EE^w+ E T.pm+p{q). {&) 

i=l i^k=l i^k=i 1=1 fc=4 i^k=\ l=i 

Here the notations il and ilm denote sums of gluon momenta ki + ki and ki + ki + km respectively. 
Expression (8) can further be simplified if we take into account that in (4) the integration over all 
momenta is done for a function which is totally symmetric in the gluons (123) and (456). Moreover, 
it is symmetric under the interchange (123) (456). Therefore, on the rhs of (8) the terms inside 
each sum are identical, and we get 

= 6P(1) - 6P(12) - 9P(14) + 9P(124) + P{q). (9) 

The 7* — rjc impact factor is given by [2] 

$^(1,2,3) = 6ei,-^0(l,2,3), i,i = l,2. (10) 

Here 

•^(l' 2, 3) = ^ M2 + (g - 2hf - ^''^ 

and = + 4m^, where is the photon virtuality and rric the charmed quark mass. The 
coefficient h is given by 

h = —ecgl-mri^hQ , (12) 
vr z 

where Cc = (2/3) e is the electric charge of the charmed quark m^^ is the r/c meson mass and 6o can 
be determined from the known radiative width T{r}c 77)= 7 KeV: 



^ = ff,/^. (13) 

3e^ V "^^c 
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The impact factor (10) is symmetric in the three gluon momenta, and it vanishes if any of the 
momenta goes to zero. Taking in (4) the product of the two impact factors and summing over 
polarizations we obtain 

F^.^^>(l,2,3)(/)(4,5,6), (14) 



where 
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^r^vc = -2 • (15) 



The pomeron P{k) can be presented as a convolution of the BFKL Green function with the 
colour distribution p{r) in the hadronic target 

Pik) = -glj d\'G{Y,k,r')p{r'). (16) 

Note the minus sign. After transforming the initial momentum space expression for the impact 
factor into the coordinate space, the impact factor is proportional to 1— exp(i/cr). When multiplying 
with the Pomeron Green's function and doing the A;-integral, there is no contribution from the '1' 
(since, in coordinate space, the Pomeron Green's function vanishes when both arguments coincide), 
and the nonzero contribution comes from the second term, — exp{ikr). The Green function has to 
be taken in a mixed representation, momentum k at the odderon side, coordinate r' at the proton 
side. Also, one side of the Green's function is amputated, the other not. 

e(^,,,,,.__L,,/_|_,.»M,(^)^" (17) 

with 

u;{u,n) = 2as{i^{l)-Rei;{^-^ + iu)), as = ^. (18) 

At large rapidities small zv's dominate. Due to the finite dimension of the target, in (17 the values 
of r' are limited by the radius R. So we can neglect factor the (r')^*'', and the integration over 
r' will be replaced by the average transverse dimension of the target R. So at high Y the Green 
function gives a factor 

\kR,^e-yeJ-'43)- (19) 



Here A = a;(0,0) = (5^4 In 2 is the BFKL intercept with a = 7asC{3). The second exponential 
factor cuts the integration over k to values In^ k < oY . However, since the integration in (4) is in 
fact convergent, we may drop this factor at high enough Y. Putting (19) into (8) we get for the 
cross-section 

i£^^,NlFPB.-\fZ,(4,\, (20) 



dt ^^"27r2g2M3 ^ ^ \ aY \M 

where the dimensionless function I{q/M) is given by the integral 

^{m)^^'I ^^^i^^'^^^)'f>i^'^^^){^\ki\-Q\Q-ki\-9\ki-k4\+9\q-ki-ki\+q). (21) 

The cross section (20) is of the order as{asNc)^ . The 8-dimensional integral (21) is non-factorizable 
and can be done only numerically. 
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As an alternative way of evaluating this contribution to the cross section, one might try to 
factorizc the integration over gluonic momenta and to transform the pomeron amplitude to the 
coordinate space using 

^^+- = -(l + 4.^)/^(^)"e^^^ (22) 
Taking the limit v ^ and introducing a function of r 




we find the integral / as 




{(6/i(r)/i(0)-6e*«^/i(-r)/i(0)-9/t(r)/i(-r)+9e^«''/i2(-r)+e^«''/i2(0)) - (r = o)}. (24) 

In obtaining this expression we used the fact that the change q ^ —q is equivalent to changing 
r — > — r in h{r). Also, we have taken into account that P{k = 0) = 0, in order to subtract the value 
of the brackets at r = and thus to improve the convergence of the integral at r = 0. Passing to 
the function 

hi{r) = h{r) - e*«''/i(-r), /ii(0) = 0, (25) 
we can rewrite (23) in a simpler form, which also shows that the right-hand side is real: 

^(f) ^(6M0)Re/.i(r)-^|^i(r)p-(l-cos(gr))/.2(0)). (26) 

Now we have only 6 integrations to be done numerically. However the presence of the oscillat- 
ing factor makes such calculations very difficult. We therefore use Eq. (21) and performed the 
integrations by Monte-Carlo methods. The results will be presented in Sec.6, together with the 
contribution from the POO vertex. 



4 The POO vertex in leading order in N(. 

Next let us consider the irreducible part Dj of the amplitude for 6 reggeized gluons. Its contribution 

to the cross section will be denoted by the superscript POO. Starting from [8] (eq.(6.3), we note 
that the rhs satisfies a BFKL-like equation which we write in the symbolic form 

{He-E)Di = DfK (27) 

Here Hq is the total Hamiltonian for 6 reggeized gluons, which is a sum of pairwise interactions 
and of gluon trajectories and describes their evolution, without changing the number of gluons. 
The energy E = 1 — j = —lu is just one minus the intercept. The driving term of the equation is a 
sum of terms which describe transitions with a change of the number of gluons, from "irreducible" 
configurations of 2, 4, or 5 gluons to 6 gluons. At this moment it is important to invoke the 
approximation of large number of colours N(. ^ oo. In this approximation, any interaction inside 
the outgoing six gluon state which connects colourless groups of gluons is damped by and can 

be neglected. This means, in particular, that once a pair of states with colour color structure of 
two odderons is formed in the driving term, Hq in Eq. (27) contains no further interaction between 
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these two odderon states. All what Hq does is to build up the bound states of the gluons (123) and 
(456). So in order to find the terms relevant for the POO transitions we only have to see whether 
the final two Odderon states couples to the driving term. One immediately sees that, at large N(., 
the irreducible configurations of 4 and 5 gluons, Dl and D^, reduce to the splitting of the initial 
pomeron into two pomerons and thus cannot couple to the two odderon final state. Therefore, the 
transitions of interest can only occur in terms which describe transitions of 2 gluons to 6 gluons. 
In [8] 4 such terms of different colour structure were found. 
The first group is given by a sum 

da,a2asda,a,aeW{l,2,3\4,5,6) +da,a,a,dasa,aeW{l,2,A\3,5,6) + (28) 

where the sum extends over all (ten) partitions of the six gluons into two groups containing three 
gluons each. Projecting onto the two odderon colour state we find the colour factor for the first 
term in (28) 

d d d d - = iVM^^°) ~ AT^ f29) 

"'aia203"'04aBa6"'ai02a3"'a4a6a6 — j\t2 — ^^c-'^c -''c > K^"^ ) 

c 

whereas for all the rest terms we have 

rl H H H _M^1)M^ /V4 .on^ 

"'aia2a4"'a3a5«6"'aici2a3"a4a5a6 — jy^2 -'^c ■ V^*-*/ 

So at large Nc wc will retain only the first term in the sum, (28). 

Apart from the W terms, the remaining driving terms in (6.3) of (27) with transitions from 2 
to 6 gluons contain 3 more groups of terms of different colour structure. Terms denoted by L in [8] 
are given by a sum 

faia2a3fa4a5ae^{^j 2, 3|4, 5, 6) + faia2a4fa3a5aeL{l, 2, 4|3, 5, 6) + (31) 

with the sum, again, extending over all (ten) partitions of the six gluons into two groups containing 
three gluons each, and the function being described in [8]. Obviously these terms give zero when 
projeted onto the colour state of two odderons. 

Finally, the terms denoted by / and J are given by sums 

d"i"~(5„5„g/(l,2,3,4|5,6) +d'^i™5(5„,„e^(l,2,3,5|4,6) + ... (32) 

and 

rf"^"^"^"Ma5a6/(l, 2, 3, 4|5, 6) + d'^^^'^-'^'^Sa^a.Iih 2, 3, 5|4, 6) + (33) 

with the sum extending over all partitions of 6 gluons into two groups with 4 and 2 gluons,known 
function / and J and 

^aia2a3a4 ^ ^^^^ai^a2^a3^a4^ ^ ^j.^^a4^a3^a2^ai^_ ^34^ 

Projecting onto the colour state of two odderons we find non-zero colour factors of the type 



d''''''''"''Sasmda,a2a3da4a,ae = 4)^ ^ ^4 _ ^^^^ 



So although these terms seem to involve transitions POO, they are down by a factor as compared 
to (29). 
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So in the end we find that, in the large Nc Umit, the transitions POO are fully described by the 
function W{1, 2, 3|4, 5, 6) which represents a convolution of the pomeron with a POO vertex. The 
functional form of W in [8] is rather complicated, however a closer inspection shows a surprisingly 
simple structure [3]. Let us briefly recapitulate this structure. It is convenient to introduce two 
operators which transform a function of momenta of two gluons into a new function which depends 
upon momenta of three gluons. Namely, define 5" to be an operator acting on 2-gluon states and 
with values on the 3-gluon states, which performs an antisimmetrization in the 2 incoming gluons, 
splits the first of them in two and sums over the cyclic permutations of the outgoing gluons: 

S{1, 2, 3|1', 2')0(1', 2') = ^ E ['^(12, 3) - </'(23, 1)] . (36) 

(123) 

Next define another operator P which performs an antisimmetrization in the 2 incoming gluons 
and splits the first of them in three outgoing gluons: 

P(l, 2, 3|1', 2'))0(1', 2') = ^[</.(123, 0)) - 0(0, 123)] . (37) 

Apart from these operators we introduce a function /(1,2|3, 4), antisymmetric in the first and 
second pairs of gluons and symmetric under the interchange (12) (34), as a sum of functions 
G(l,2,3), which were introduced in [9, 10] in the context of the three-pomeron vertex: 

/(1, 2|3, 4) = G(l, 23, 4) - G(2, 13, 4) - G(l, 24, 3) + G(2, 14, 3) . (38) 

The explicit form of the general function G(l,2,3) is not important for our purpose (it can be 
found e.g. in [9, 10]). We only have to know that 

G(l,2,3) = G(3,2,l), G(0,2,3) = G(l,2,0) = 0, (39) 

and that, up to a coefficient, G(1,0, 3) is given by the BFKL Hamiltonian H2 applied to the 
pomeron: 

G(l,0,3) = -i-(i?2P)(l,3). (40) 
In terms of S, P and / we find 

W{l,2,3\4,5,6) = -pt{Si-Pi)fu{Sl-P^), (41) 
where the indices 1 and 2 refer to the triplets of gluons (123) and (456), respectively. 



5 Part of the cross-section with a POO transition 

To find the cross-section corresponding to the rjc production via the POO transition (Fig. 1) we 
have to couple the POO vertex with the two odderons attached to the initial and final 7* rjc 
impact factors. To write it in a compact form we introduce the Green functions Gg^^ and Gg^^ for 
the initial and final odderons composed of gluons 123 and 456 respectively. They evolve the odderon 
state in the rapidity interval with length y <Y. Then using (41) we can write the cross-section as 



dt 



= -LjN^Fi^oo)^ [dy< ^\\gS\Si - Pi)fu{Sl - P^)G^^^M > , (42) 
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where it is assumed that the averaging is done independently over the gluons 123 (index 1) and 
456 (index 2), function / depending on both groups of variables. 

At this point we recall that the full set of odderon states and consequently the full Green 
function G3 are unknown. We are going to use a part of it corresponding to the solutions found 
in [3], which have the maximal intercept of all known states and besides have a non-zero coupling 
to the perturbative 7* rjc impact factor. These odderon states are expressed via the known 
antisymmetric pomeron states E^'^'"'^ with odd values of n: 

^'('''-)(1, 2, 3) = ciu, n)-^S{l, 2, 3\l' ,2')k[\''E('^'-\l' , 2') , (43) 
where 



c(i,, n) = J — rA?^r T (44) 

^ ' ' Y -3(27r)3cj(zy,n) ^ ' 

and io{i',n) is given by Eq. (18). The part of the Green function G3 corresponding to these states 
then aquires a form similar to the pomeron Green function [5] 

G3(2/i|l,2,3|l',2',3') = J2 yt^i^e^'''^'''"^/3(j^,n)*^'''''ni>2,3)*(^'")*(l',2',3'), (45) 

odd n 

with 

- (27r)2(;/2 + nV4) 

~ [1/2 + (n-l)V4][z^2 + (^+ 1)2/4] (46) 

and ^'(''''*) given by (43). 

Now we use the fact that we only study our cross-section in the region where both the rapidity 
of the pomeron y and that of the odderon yi = Y — y are large. This allows to retain in (45) only 
the branch |n| = 1 with a maximal intercept and also restrict the integration over v to small values. 
In the limit u ^ Q the coupling of the odderon to 7* r)c impact factor was calculated in [5] to 
give 

(^^^.J*^'''")) = --&%-^T^'^'"^Xa72 • (47) 
' ' TT q c[i>, n) q'^ + 

After summation over the photon polarizations the two matrix elements (47) provide a factor 

1 ^ ,q^'^'^^-^-\ (48) 



TT^ c(z/i,ni)c(f2,n2) (g2 -|- M^)^ 

where {vi,ni) and {^2,^2) refer to the summation and integration variables in and ^3'^^ 

respectively. At finite value of q one can neglect the last factor in (48). 
We are left with the matrix element 

T ^< *^"^)|(5i - Pi)/i2(5| - P2^)|$^'"^) > . (49) 

Its calculation is greatly simplified by the property of the odderon state (43) found in [3] . For any 
function 0(1,2) of two gluon momenta and odderon state (43) one has 

< ^M^^g _ p)|^ _J_ ^ ^(</,n)|^ > _ (50) 



12 



Note that the matrix element on the left-hand side is taken in the space of three gluons, whereas 
that on the right-hand side is taken in the space of only two gluon momenta. This property greatly 
simplifies the matrix element (49). Using (50) we get for it 

T = ^- < S^'"^Vi2|£^^'"^^ > ■ (51) 

Now we again use the fact that yi is large and so only values |n| = 1 and \v\ « 1 contribute. 
At \n\ = 1 and small v the pomeron wave funtions entering (51) reduce to 5 functions of gluon 
momenta [5] 

E^^'^'-\l,2),^o = ^('^'(^i) - <^'(^2)), £;('^'±^)(4,3),_o = ^('^'(^4) - 5\k^)) . (52) 

Note that the second wave function has to be taken conjugate. Putting this into (51) and recalling 
(38) and the properties of the function G (Eq. (39) and (40) we get 

T=—. ^, -—^G{q,Q,-q) = A^±(H2P){q) . (53) 

c{yi,ni)c{v2,n2) (27r)^(?^ c(i^i, ni)c(i/2, n2) tt^^ iVc 

The last factor is just the BFKL Hamiltonian applied to the Pomeron state. 

As in Sec. 2, to find the pomeron P{q) attached to the hadronic target we present it as the 
BFKL Green function applied to the colour distribution in the hadron, Eq. (16). Again we need a 
mixed amputated-non-amputatcd Green function in the momentum space in its amputated part, 
Eq. (17). Applying the BFKL Hamiltonian we get 

H2Giy, q, r) = ^qr J duu{u, 0)e^-(''''^) . (54) 

At large y with finite q and r we neglect all dependence on u except in the exponential to obtain 
similarly to (19) 

HG{y, q, r) = J^grAe^^. . (55) 
27r^ y ay 

Integration of r with the target colour density converts it into the average target transverse radius 
R with a minus sign , see (16)). 

So we find for the matrix element T 

T = 1- ELqR/^ey^ JT, (56) 

27r^g^ c(z^i,ni)c(z/2,n2) A'c V ay 

We have finally to do the integrations over vi and V2 and summations over ni and n2, which in the 
limit of high yi only take values ±1. These latter summations are trivial and give a factor 4. At 
|n| = 1 and small v we have 

a;(z/,±l) = -2a,C(3)z^^ (57) 

One of the denominators in (46) reduces to and is singular at i/ — >^ 0. However this singularity is 
cancelled by the square of the \l(?(y^ ±1) coming from (48) and (56). Therefore we find at small v 



C2(l/,±1) 



= 127r^C(3) . (58) 
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Neglecting all the rest i/-dependence except in the exponential in (45), integrations over vi and 1^2 
provide a factor 

(59) 



2a.C(3)yi ■ 

Combining all the factors we finally get for the cross-section a simple expression 
at J Ncq(q^ + M^Y agViV ay 



Ncq{q^ + M'^Y asViM ay ' 

It steadily grows as diminishes and behaves as 1/g at small q. Integrating over all q we find the 
cross-section 

^ " " J 'N, M3 asViy ay ^ ' 

It falls with and the meson mass as 1/(Q^ + mp5')3/^. The cross-section has an order as{cisNcf' , 
an order higher in agNc than the leading contribution given by the reggeizing part (Fig. 2). This 
implies that the vertex POO has the same order agNc as the triple pomeron vertex. 



6 Numerical results 

Both the cross-section with a pure pomeron exchange and with a POO vertex have a simple de- 
pendence on energies, which separates into a factor with the expected behaviour at large rapidities. 
Separating also all the rest non-trivial factors we find the pure pomeron exchange contribution as 

^ = c(-)ae.a.(c..)^6g^/ (£) f^^\Y) , (62) 
r\Y)=e'^\[Z (63) 



dt 

where 
and 



c(^) = (64) 

5184 7r6 • 



The part originating from the POO vertex is 



= c^^oo^a^mas J dy i-sM ^JI^^^^J '^^^Hy, y) , (65) 
where now 

/^^°^H>^,y) = -7T^e^^./^ (66) 
as{Y-y) y ay 

and 

^(POO) _ ^i^^^^C(3)ln2 
The cross-section with the POO vertex integrated over all transferred momenta is 

^(POO) ^ l^e(^oo)^^^^^ J (a,)662!^/(^oo)(y, y) . (68) 
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With TVc = 3 the colour factors become 

fP = 200 • 3-^, = 1600 • 3-^ . (69) 

Regarding the region of integration in rapidity y, as explained in the introduction, we consider the 
interval 6Y < y < Y — 6Y to warrant the use of the asymptotic forms for both the pomeron and 
the odderons. We choose SY = 3. 

One should also take a certain care with the coupling constants in the expressions for the cross- 
sections. In fact they refer to different scales relevant to the studied processes. Obviously one of the 
coupling constants refers to the coupling to the proton at a small and so non-perturbative scale. For 
the process mediated by the pure pomeron all other coupling constants are to be taken at the scale 
of the 7* — ^ rjc transition, that is the maximal of q and M. The cross-section falls quite rapidly as 
q becomes larger than M, so that we can safely take M as the relevant scale. For the process with 
the POO transition however only three of the remaining 6 a's refer to this scale. Other three are 
to be taken at an intermediate scale, characteristic for the POO transition at rapidity y. For high 
rapidities Y and y one can use the fact that the characteristic momenta k in the BFKL pomeron at 
rapidity Y have the order Ink ^ VY. Then one obtains a crude estimate for the coupling constant 
at the POO junction as 




apoo ^ \l —(^siM) . (70) 
We have taken Y = ln(l/a;) with x defined as 

Passing to concrete values of the coupling constants we take as{M) as given by the leading 
order /3-function with 3 or 4 flavours Nf and Aqcd = 0.2 GeV/c. The value of ccg at the POO 
junction was taken according to (70). As to the values of the pomeron intercept and its coupling 
to the proton, we have borrowed them from [11], where the proton structure function at small x 
was fitted by the pomeron exchange. From this fit one extracts both A and the product agR: 

A = 0.377, a^i? = 0.096 fm, {Nf = 3), a^i? = 0.058 fm {Nf = 4) . (72) 

At first sight one may expect a large difference in the results for different Nf. However a smaller 
value of agR for Nf = 4 is compensated by a larger value of the rest of the coupling constants, so 
that the final results are practically independent of the number of flavours taken into account (see 
Fig. 7 below). 

The calculation of the POO contribution (65) is straightforward. To find the contribution from 
the pure pomeron exchange (62) one has to calculate integral (21). We did this using the standard 
Monte-Carlo program VEGAS. The results for I(q/M) are presented in Fig. 4. 

Our final results for the differential cross-sections (62) and (65) and their sum (for Nf = 3) are 
presented in Figs. 5 and 6 for ^/s = 300 GeV and = and 25 (GeV/c)^ respectively. 

As we see, the contribution from the POO transition turns out to be of the same order as the 
one coming from the direct coupling of the pomeron to non-interacting gluons (P contribution). 
However the two contributions seem to behave differently at small transferred momenta. From (65) 
we see that at g ^ the POO contribution rises as 1/q. On the other hand, the P contribution 
does not show such a behaviour and seems to tend to a constant or zero at small q. This may help 
to see the POO contribution against the reggeizing pomeron contribution at very low q. 
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Figure 4: The function I[q/M) from eq.(21). 
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Figure 5: Differential cross-sections da/dt from the reggeized pomeron exchange (P), POO transi- 
tion (POO) and total at Q"^ = 0. 
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Figure 6: Differential cross-sections da/dt from the reggeized pomeron exchange (P), POO transi- 
tion (POO) and total at =25 (GeV/c)^ 



Integrated over transferred momenta the total cross-sections are found at = to be 
a(^) = 34pb, a(^00) = 31pb, a = cr^^) + a^^^O) = 65 pb 
and at =25 (GeV/c)^ 

a(^)=0.87pb, a(^OO) = 0.67 pb, a = a^^) + a^^^^) = 1.54 pb . 

The total integrated cross sections a (sum of P and POO) in the range of photon virtualities 
< < 25 GeV~^ are shown in Fig. 7. Here we presented results for both Nf = 3 and 4. As one 
observes the difference is quite insignificant. 

With the growth of energy both contributions increase, preserving their shape in q. The increase 
is much more pronounced in the reggeizing pomeron contribution, since in this case the pomeron 
occupies the whole rapidity range, whereas for the POO transition this range is shorter. 

7 Conclusions 

We have calculated the cross section of inclusive diffr active photo- and leptoproduction of r]c mesons, 
7*P Vc+^- We have considered the 'triple Regge' contribution which contains the coupling POO 
of the pomeron to two odderons. Inclusion of a second contribution where the pomeron directly 
couples to two three-gluon states results in a significant rise of the cross-section which grows with 
energy. However, in order to see the structure of the QCD odderon state with C = — 1 state one 
has to select diffractive events with a large enough gap between the missing mass state "X" and 
the r]c- The total production rate is found to be of the order 60 pb for photoproduction. 

In our previous publication [5] we calculated the production rate for the quasielastic reaction 
^y* + p ^ rjc + p due to odderon exchange. The photoproduction cross-section was reported to 
be 27 pb, which would be smaller by a factor of 2 compared to the present case. However, note 
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Figure 7: Integrated cross-sections a for different number for flavors as a function of Q^. 



that for the quasielastic process considered in [5] one had to make an assumption of the non- 
pertTirbativc oddcron-proton couphng. In [5] wc used the coupUng proposed in [2], and we put the 
effective couphng constant equal to unity. However, in a recent analysis of the pp and pp elastic 
scattering data this coupling constant has been estimated to be 0.3 [12]. With this value of the 
effective coupling constant the cross-sections reported in [5] have to be reduced by a factor 30, and, 
in fact, the cross section of the inclusive rjc is much larger than the quasielastic one. 

In the present calculation the very poorly known odderon-proton coupling does not enter. In- 
stead one has to know also the non-perturbative pomeron-proton coupling, transformed into the 
value of the product a^i? where R is the effective proton radius. This product can be found with 
a much higher degree of reliability. In this study wc have used the fit to the experimental proton 
structure function in [11]. Note that it gives physically reasonable values A = 0.377 and R = 0.59 
fm (for Nf = 4), which more or less agree with estimates made by different methods. Correspond- 
ingly, we feel that the cross-sections found in this paper are much less affected by the uncertainty 
in the nonperturbative coupling of the proton. 
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